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Abstract 

This paper addresses the existence and uniqueness of strong solu- 
tions to the stochastic variational inequality 

m 

dX — AX dt + F(t, £, X)dt + fi(X)dt 3 ^ X ^ k d^ k {t) + / \t)dt 

k=l 

in a bounded domain O cM. d with Dirichlet homogeneous conditions. 
Here /3(r) = for r > 0, /3(0) =] - oo, 0], 0(r) ^ for r < 0. An appli- 
cation to the one-phase Stefan problem with a stochastic heat source 

is given. One studies also the corresponding stochastic parabolic equa- 

dX 

tion with Signorini boundary conditions — h 0(X) 3 on dO. 

ov 
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1 Introduction 

This work is concerned with the existence and uniqueness of the stochastic 
variational inequality 
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' d t X(t, - \X(t, Odt + F(t, e, X(t, 0)dt 



rn 



< 



in (0,T) x = Q T , 
^ X(t,£) = on (0,T) x = £ T , 



fc=i 



(1.1) 



where (9 is an open and bounded domain of M. d , d > 1, with smooth boundary 
9(9, {^(t)}^! is a system of mutually independent Brownian motions on a 
probability space {Q, J 7 , P} and /3 : R — > 2 R is the multivalued function 

(3(r) = for r > 0, 0(0) =] - oo, 0], f3(r) = for r < 0. (1.2) 

As regard the functions F : Q T x H. — >■ R, ^a, : [0,T] x -> 1, we assume 
that the following hypotheses hold. 

(i) F = F(t,£,r) is measurable in (t,g) G Q T , F(t,£,0) = and 



\F(t,£,r)-F(t,Z,r)\<a\r-r\, Vr, f el, a.e. (U) £ Qt (1.3) 



The function / : Qr x f2 — )■ R is an adapted process from (0, T) to L 2 (0) 
with respect to the probability basis J 7 , {J-^}, P}, which will be made 
precise later on. 

Formally, problem (1.1) can be rewritten as a free boundary value problem 
(the stochastic reflection problem) 



where a > 0. 



(ii) fi k e C 2 



([0,T] x O), k = l,...,m. 
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d t X - A ( X dt + F(t, f , X)dt = / dt + ^ X/j fc (i/3 fc 



in [(U)eQ r ; *(*,£)> 0] 



(1.4) 



< 



X > 0, d t X - A^X dt + F(t, f , X)dt > fdt + Y^ Xii k dp k in Q T 



k=i 



X(0) = xinO, X = on E T . 
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As regards the literature on stochastic parabolic variational inequalities of 
this form or, more generally, on parabolic differential equations driven by 
Gaussian noise and reflected to the boundary of a convex set of R d , the 
works [6], [8], [10], [11], [13] must be cited. (See also [2], [3], [7] for the finite- 
dimensional case and [4], [5], [14] or [1], Chapter 1, for existence theory in 
infinite-dimensional Hilbert spaces.) It should be said, however, that all refer 
to the case of stochastic equations with additive cylindrical Gaussian noise 
and the solution X is taken in a generalized sense. Roughly speaking, such a 
solution X satisfies (1.1) where (3(X) is replaced by a measure on (0, T) x O 
with the support in [X = 0]. Theorem 2.2 below provides the existence and 
uniqueness of a strong solution to (1.1) which is more appropriate of classical 
formulation (1.4) of the obstacle parabolic variational inequality. We develop 
here a parallel study for the stochastic parabolic equation (1.1) with Signorini 
boundary value conditions. Namely, 



f d t X(t, - A ( X(t, Odt + F(t, £, X(t, 0)dt 

m 

= f{t,Odt + J2x(t,OVkd{3 k {t) inQ T 



k=i 



3X 



(1.5) 







where — — is the outward normal derivative to dO. 
av 

As made precise in Section 4, problem (1.1) can be used to describe the 
dynamic of the one-phase Stefan melting (solidification) problem in presence 

m 

of a stochastic heat source of the form ^^#/i fe <i/? fc , where 9 = 0(t,£) is the 

k=l 

temperature. As regards (1.5), it is a stochastic parabolic equations with 
unilateral conditions on the boundary. 



2 The main results 
Notation 

Given the stochastic basis {Q, P, J 7 , {J^j^o} and Z a Hilbert space with the 
norm ||-|| z , we denote by Mp(0, T; Z) the space of all progressively measurable 
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Z-valued processes X : Vt x (0, T) — > Z with respect to the filtration {F t }t>o 
such that 

where E is the expectation. Denote by Cp([0, T]; Z) the space of all processes 
X G C([0,T];L 2 (n,F,J,Z)). By L p ad (Q; C([0, T]; L 2 (C)) we denote the 
space of all adapted processes with respect to {J 7 *}, X : — >■ C([0, T]; L 2 (C)) 
such that 

^ll^llc'([0,T];L 2 (Cl)) < °°- 

By L p (0), where (9 is a bounded and open subset of R d and 1 < p < oo, we 
denote the space of all p-summable functions on O with the norm denoted by 
| ■ | p . By L p (0, T; Z) we denote the space of all Bochner Z- valued p-summable 
functions on (0, T) and by W 1,p ([0,T]; Z) the space of all absolutely conti- 
nuous functions y : [0, T\ — > Z such that — ^ G L p (0,T; Z). The norm of Z is 
denoted by || ■ \\z- 

In the following, Hq(0) and H 2 (0) are standard Sobolev spaces 

on (9. We also denote by H~ l (0) the dual space of Hq(0) and set 

# 2,1 (Qt) = \y e L 2 (o, r ; 1^(0) n H 2 (0), ^ g L 2 (Q T )j . 

Definition 2.1 By strong solution X to (1.1), we mean a process X G 
Cp([0, T]; L 2 (C)) n M|(0, T; ^(C)) such that there is 77 G M|(0, T; L 2 (C)) 
and P-a.s., 

' ^(U)eWU)) a.e. (U)eQr 

X(t)=x + [ A^X(x)ds-[ F(s,£,X(s))+ri(s))ds 
< Jo Jo (2.1) 

+ / f(s)ds+y2 / A*fc(s)^(s)d/3 fc (s), tG[0,T], 
I J ° k=i Jo 

where the integral is taken in the sense of Ito in the space 11^(0). (We refer 
to [9] for existence and uniqueness of such a solution.) 

It should be noted that, for any such a solution X to (1.1) in the sense of 
Definition 2.1, the Ito's formula is applicable. 
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Taking into account that the relation r) G B(X) is equivalent with 

X > 0, r/ > 0, Xr) = 0, a.e. in Q T , (2.2) 

it follows by (2.1) that X defined above is a solution to (1.1) in the sense 
of (1.2). 

Theorem 2.2 Let f G M|(0, T); L 2 (0)) and x G H^{0) be such that x>0, 
a.e. in O . Then (1.1) has a unique strong solution X . Moreover, we have 

X G If ad (n- C([0, T}; L\0)) n M£(0, T; if 2 (G)), (2.3) 

/or all 1 < p < 2. 

In particular, it follows by (2.3) that A 5 X G L 2 (Q T ), P-a.s., and so (2.1) 
holds, a.e. on (0, T) x O x ft. 

Now, we come back to (1.5) to define a strong solution for this equation. 
To this purpose, we define the multivalued nonlinear operator A(t) : V — > V, 
V = H\0), 

v .(A(t)X,<p) v = f(VX(Z).V<p(t) + F(t,t,X(Z))<p(Z))dt 

J ° , (2.4) 
+ / r,(X)(Z)<p(Z)dZ, V(/3 g V, te[0,T], 

JdO 

where r\ G L 2 (dO) and 

TjpOtf) G P(X(0), a.e. (2.5) 

Here V is the dual of V and y/ (•, -) y is the duality between V and V induced 
by the space H = L 2 (0). (We have V C H C V in algebraic and topological 
sense.) 

The operator ^4 is not everywhere defined on V but on the set D(A) of 
all X E V for which there is rj G L 2 (dO) satisfying (2.5). Roughly speaking, 

A(t) = —A with the nonlinear boundary conditions — h B(X) 3 on dO. 

ov 

Definition 2.3 X is said to be strong solution to (1.5) if the following con- 
ditions hold: 

X G C P ([0, T\- L 2 (0)) n M P 2 (0, T; V), (2.6) 
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X(t)-xe - f A{s)X{s)ds+ [ f(s)ds 

J ° J ° (2.7) 



-a.s. 



+ / X(s)J2^(s)dp k (s), Vte [0,T], E 

J ° fc=l 

Here / A(s)X(s)ds is the Bochner K'-valued integral of the multivalued 
Jo 

function t A(t)X(t) and / X(s) V fi k (s)d(3 k (s) e L 2 (0), P-a.s., is taken 

Jo k=i 

in the Ito sense. 

We may rewrite (2.7) as the multivalued stochastic integral equation in V 

m 

dX(t) + A{t)X{t)dt 3 f(t)+J2 X (t)vMt), t G (0,T), 

k=i ( 2 - 8 ) 

X(0) = x. 

Theorem 2.4 Let f e M$(0, T); L 2 (0)) and x e H%(0) be such that x>0, 
a.e. in O. Then (1.5) has a unique strong solution X. Moreover, X > 
a.e. in (0,T) x O x fl 

In particular, it follows by (2.3) that (2.7) holds, a.e. on (0,T) x O x fi. 

3 Proofs 

3.1 Proof of Theorem 2.2 

We set 

m 

and 

fc=i ^ ' 
By the substitution 

X(t,0 = e^y(t,0, cuen, (3.1) 
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equation (1.1) reduces to the random parabolic equation 

f ^ (t, - e-^A(e^y(t, 0) + e^F(t, £, e^y(t, 0) 

+Jl(t, £)y(t, + /3(y(t, 0) 3 /(*, 0, (*, e <5t, (3.2) 
2/(0,0 = ^(0 m 0, 
[ y(t,0 = on E T . 

Indeed, if y = y(t,0 is a sufficiently regular solution to equation (3.2), then, 
by Ito's formula, we have 

dX = yd(e^) + ^ dt in Q T 

and 

N 



d(e") = e»dn+le»J2fi 2 k dt 



2 

k=i 

N N N„ 

= e*Y,to d fc+ 2^12 titi + e*Y, -Iff fa dt - 

k=l k=l k=l 

Then, substituting X into (1.1), we obtain (3.2), as claimed. 
Now, by (3.2), we have P-a.s. 

^-Ay + F(t,y)+g-Vy + /3(y)3f in Q T , 

y(0) = x in O, 
y = on Ey, 

F(t, y) = e-"(«)F(t, £, e**.%) + /ly - (| V^| 2 + A/i)y, 
y = -2V/x. 



(3.3) 



where 



(3.4) 



Any regular and predictible solution t — > y(t) with respect to the stochas- 
tic probability basis {Q, J 7 , J~t, IP} to equation (3.3) leads via transformation 
(3.1) to a solution X to (1.1) in the sense of the above definition (see, e.g., 
[9], p. 232). 
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We note that (3.3) is a random parabolic variational inequality for which 
the standard existence theory applies (see, e.g., [1], p. 209). Namely, for each 
uj G f2, (3.3) has a unique solution y = y(t,£,uj), 

y G W 1 ' 2 ([0,T};L 2 (O))DL 2 (0,T;H^(O)r)H 2 (O)), (3.5) 

which is obtained as limit for e — > of the solution y e to the penalized 
equation 

dy ~ 

-^-Ay £ + F(t,y £ )+g-Vy £ + f3 £ (y £ ) = f in Q T , 
y £ (0) = x in O, y £ = on £ T , 

where (3 £ (r) = - (r - (1 + £/3) _1 r) = -- r~, Vr G R. 

Indeed, if y £ G ^([O, T]; L 2 (C)) n L 2 (0, T; ^(C) n H 2 (0)) is the solu- 
tion to (3.6) taking into account that, by (1.3), a.e., (t, £) G Qt, 

|F(t,r)| <a|r|, Vr G R, (U) 6 Qt, 

|s| < Csup{|&(t)|, te[0,T], fc = l,..,m} = ^), wGfi, 
after some straightforward calculation, we obtain the estimate 

\VMI + f U(s)\\ 2 w ds < \x\ 2 2 + jf \f{s)\\ds + 5 2 (u), t G [0, T]. (3.8) 

Similarly, multiplying (3.6) by f3 £ (y e ) and integrating on (0,t) x (9, we 
obtain 



J\\(3 £ (y £ (s))\l + \Ay £ (s)\ 2 2 )ds < C Qf \f(s)\ 2 2 ds + t?( u )) , 



(3.9) 
WG [0,T], 



because 



/ c ¥ (W ^ = ^ 1 a ' 6 ' tG [0 ' T] ' anda " = in °°- 

(Everywhere in the following, we denote by C several constants independent 
of e and G fi.) Taking into account that E(5 2 ) < oo, we obtain by (3.8), 
(3.9) that 
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E 



T 



lv*(t)l! + / (\\vMfkw 



+ 



dye 
dt 



(*) 



+ \\Ay £ (t)\\t + \P £ (y £ (t))\ 2 2))dt 



(3.10) 



<c(\x\l + E£\f(t)\ 2 2 d?j, Vte[0,T}. 



We note also that, since y £ is obtained by an iteration process given by 
the contraction principle, y £ (t,u) is progressively measurable and, therefore, 
adapted to the filtration {J-t}t>o- 

Moreover, by (3.6) we have also for e,e' > 



1 f* 

-\y E (t)-y E ,(t)\ 2 2 + / \V(y £ (s)-y £ ,(s)\ 2 ds 



+ 



[ [ My £ (s)) - p £ ,(y £ (s)))y £ (s) - yAs))dsdt 
Jo Jo 

<\f Jjdw g\\y £ (s)-y £ ,(s))\ 2 dsdt 

<sup{|div#|; / \y £ (s) - y £ ,(s)\ 2 2 ds 

Jo 



and 



and so, by (3.8), (3.9), we have 

\Ve{t)-yAt)\l+ f \V(y £ (s)-y £/ (s)\ 2 2 ds < C^is + s'), Vt e [0,T\. (3.11) 
Jo 



By (3.10), (3.11), it follows that there is an adapted process 
y e M P 2 (0, T; H^O) H L 2 ad (Q; C([0, T]; L 2 (0)) 
^M P 2 (0,r);I 2 (O)), 



(3.12) 
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(3.13) 



such that, for e — > 0, 

y £ -+ y strongly in L 2 ad (n;C([0,T]; L 2 (0)) H M 2 (0,T; H^O)), 
and P-a.s., 

Ay £ ->■ Ay weakly in L 2 (Q T ), P-a.s., 
->■ ^ weakly in L 2 (Q T ), P-a.s., 

F(i,j/ e ) -> F(*,y) strongly in L 2 (Q T ), 

^ % weakly in L 2 (Q T ), 

where r\ G (3(y), a.e. in Q, P-a.s. We have also P-a.s. 
dy ~ 

■^-Ay + F(t,y)+g'Vy + ri = f, a.e. in Q T , 

2/(0,0 = ^(0, in 0, 2/ = on S T . 

Since y is a smooth y G W /1 ' 2 ([0, T]; L 2 {0)) and predictible solution to (3.14), 
we conclude that X given by (3.1) is the strong solution to (2.1) (see [9], 
p. 217). The uniqueness of the solution X follows by a standard argument 
involving Ito's formula. Moreover, since 

E{sup(e^ + |A/x|; t G [0, T]) n } < oo, Vn G N, 

we infer by (3.12) that (2.3) holds. This completes the proof. 

3.2 Proof of Theorem 2.4 

Proceeding as in the proof of Theorem 2.2, we reduce, via transformation 
(3.1), equation (1.4) to the random equations 

r dv ~ 
-^-Ay + F(t,y)+g-Vy = f in Q T , 

< j/(0) = x in O, (3.15) 

l! + !v + /»<»)30 onE T , 

where F and y are given by (3.4). Now, we fix uj G Q and treat (3.15) as a 
deterministic nonlinear parabolic equation. 
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It must be said that, though the existence theory for nonlinear parabolic 
boundary value problems is not directly applicable in the present situation, 
we may prove, however, existence in (3.15) via standard approximation tech- 
nique. Namely, we approximate (3.15) by the penalized equation 



-^-Ay £ + F(t,y £ )+g-Vy £ = f in Q T , 

y £ (0)=x in O, 

dy £ dfi 



(3.16) 



Equivalently, 



dy 



on Tit- 



-£+Mt)Ve = f(t) in (0,T), 



y e (o) = x, 

where A £ (t) : V — > V is defined by 



(3.17) 



v 



(Vy £ ■ Vip + F(t, y £ )y + (g ■ Vy £ )ip)d£ 

(3.18) 



+ 



We have, for all t G [0, T], 

\\A(t)y\\v> < C\\\y\\ v , VyeV, (3.19) 
v ,(A e {t)y,y) v > C 2 \\y\\ 2 v - C 3 \\y\\ 2 H , Vy G V, (3.20) 

v ,(A £ (t)y-A £ (t)y,y-y) v > -C 4 \\y - y\\ 2 H , Vy G V, (3.21) 

where Q, i = 1,..,4, are independent of e. To get (3.20), (3.21), we have 
used the trace-interpolation estimate 



no 



- c ^U(do) - ^II^IU^oII^IIhho). 



Then, by a standard existence result for the Cauchy problem (3.17) (see 
J.L. Lions [12] or [1], p. 177), there is a unique y £ G L 2 (0,T; V)nC([0, T];H), 

with ^ G L 2 (0,T;V') which satisfies (3.17), a.e. on (0,T). It should be 



11 



mentioned that, since the solution y £ can be obtained via Galerkin scheme 
or by iteration via Banach contraction principle, the process y £ = y £ (t,u) is 
adapted in the stochastic probability basis {f2,P, J 7 , {J~t}t>o}- 

Now, by (3.17) or by (3.17), (3.18), we have the following a priori esti- 
mates 

\ J t be{t)\\ 2 H + l|Vy e (f)||!r < C(||y e (*)||lr + WvMW^o)) + H/WIlL 

a.e. t e (0,T). 

This yields, as in the proof of Theorem 2.2, 

E\\y £ (t)\\ 2 H + E [ \\y £ (s)\\ 2 v ds<C, t e [0,T]. (3.22) 
Jo 

(We shall denote by C several positive constants independent of e.) 
Now, we take in (3.3) ip = j3 e (y e ) and obtain that 



4 / je(y £ (t,OW+ I \fc(ye(t,m 2 dt<c(\\ y M\l + \\f(tn 2 Hl 

O JdO 



dt 



a.e. t E (0,T), 

where j £ (r) = / /3 £ (s)ds. This yields 
Jo 

[ j £ (y £ (t,OH+ f [ f3 2 £ (y £ (s,OH 

Jo Jo Jao (3>23) 

< / j e (x)d£ + C [ \\f(s)\\ 2 H ds<C,Ve>0,te[0,T\. 
Jo Jo 

Finally, if in (3.22) we take ip = y £ — y £ i and use (3.17), we obtain that 

+ f [ ((3 £ (y £ (s,0)-Ps'(ye'(s,0))(ys(s,0-ys'(s,0)dsdt 

Jo JdO 

<C(co) [ \\y £ (s)-yAs)\\ 2 H ds, Vte[0,T]. 

Jo 
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Then, by (3.9) we see that 

\\y £ (t) - y £ ,(t)\\ 2 H + [ \\y e (s) - y £ <(s)\\ 2 v ds <C(e + e')y, P-a.s. 
Jo 

and, therefore, 

y e -+ y in L 2 (0, T; V) D C([0, T]; H) (3.24) 

as e — > 0. Note, also, that the process y = y(t, ui) is adapted. Moreover, since, 
by (3.9), {(3 £ (y £ )} is bounded in £ 2 (£), we have (eventually on a subsequence) 

/3 £ (y £ ) ->■ 7/ weakly in L 2 (S T ). (3.25) 

On the other hand, by (3.24) and the trace theorem, it follows that, for 
e ->■ 0, 

y e -> y strongly in L 2 (£ T ) . (3.26) 

Since f3 is maximal monotone in L 2 (£^) x L 2 (St), we infer by (3.25) and 
(3.26) that 

a.e. (U)eEr. (3.27) 
Then, letting e tend to zero into (3.3), we see that 

! + *)» = /«), a.e. t 6(0,T), (32g) 
1/(0) = 
where is given by 

v>(A(t)y,ip) v = I (Vzy)-Vv + F(t,y)v + (g-Vy)v)dt 
Jo 



90 



(3.29) 



Now, if we take X = e^y as in (3.1), taking into account that y is a 
regular predictible solution to (3.28), we conclude that X is a strong solution 
to (1.5) in the sense of Definition 2.3. By (3.23) it follows also that 

y(t, C, u) > 0, a.e. (t, f, uj)eQ T x n, 

and, therefore, X > 0, a.e. in Qt x ^- The uniqueness follows by (2.7) 
(or (2.8)) via Ito's formula, but the details are omitted. This completes the 
proof. 
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4 The stochastic one-phase Stefan problem 

The Stefan problem describes the heat conduction in a medium O C M. d , 
d = 1,2,3, phase change. Assume that, at time t > 0, O l C O is the solid 
region (ice) and that O* = {£ e 0\ < £{£) < t < T}, dO l = {t = £(£)}, 
where £ is a smooth function. Then the classical one-phase Stefan problem 
is described by the system (see, e.g., [1], p. 231) 

' ^-A9 = in {(U)gQt; £(0<t<T}, 

/ = in {(*,£); < < (4.1) 

v i e-ve = - P on 5 = {(t,0; * = 

k f(0,0 = ^o(0 in 0, = on S T , 
where 6* = 0(t,£) is the temperature. Under a stochastic multiplicative per- 

m 

turbation of the form W = 9^^ fikfik, system (4.1) becomes 

fc=i 

m 

d9-A6dt = eJ2vkd{3 k in {(t,£) e Q T ] <t <T}, 
k=i 

9 = in {(*,£); 0<t <£(£)}, (4.2) 

V f 6» • W = -p on 5, 

0(0) = O in O, = on E. 

Roughly speaking, this means that in (melting) processes arises a heating 

m 

noise source of the form ^^9/ik/3k- 

k=i 

We shall show below that equation (4.2) reduces to an obstacle stochastic 
problem of the form (1.1). To this end, it is convenient to reduce (4.2) to a 
random one-phase Stefan problem via the transformation 

9 = ef i z in Q T . (4.3) 

Arguing as in the proof of Theorem 2.2, we see that this happens if z is 
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solution to the equation 

dz 



dt 



- Az - 2V/i • Vz + (n - Afi - | V[i\ 2 )z = 



in {(U)eQ T ; £(0<t<T}, 
2 = in {(t,0; o<t<e(t)}, 

V € 2 • W = -p on 5, 
k z(0) = 6> in O, z = on S T . 



(4.4) 



Equation (4.4) is a random one-phase free boundary parabolic equation of 
Stefan type which can be reduced to a parabolic variational inequality via 
the transformation 



/ ft 



y(t,0 = < 



[ 9(s,S)ds if ieO\0^ t>£(0, 
[ 6(s,£)ds if f G O , <t <T. 



We set 



f(t,0 = 



-p if f G < i < r, 

O (0 if 0<t<T, 



where C° = {£ G 0\ O (O > 0}. 

Then, arguing as in Lemma 5.1 in [1], it follows that y is the solution to 
the equation 

^ - Ay - 2Vp ■ Vy + (£ - - |V/i| 2 )y = / O l 0+ in P'(Q), g 
y(0) = in £>; y = on E, 

where Q + = {(£,£); < t < T} and V(Q) is the space of distributions 
on Q. Equivalent ly, 



5y 

at 



- Ay - 2Vp • Vy + (fJ. — Afj, — |V^| 2 )y + /?(y) 3 f l Q+ 



y(0) = in O; y = on E r , 



(4.6) 



and so, X = e M y is solution to equation (1.1), where / = folQ + e fl . In this 
way, the existence into problem (4.2) reduces to Theorem 2.2. 
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